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Abstract 

The asymptotic distance between trajectories g?oo, is studied in detail to char- 
acterize the occurrence of chaos. We show that this quantity is quite distinct 
and complementary to the Lyapunov exponents, and it allows for a quantitave 
estimate for the folding mechanism which keeps the motion bounded in phase 
space. We study the behaviour of cfco in simple unidimensional maps. Near 
a critical point doo has a power law dependence on the control parameter. 
Furthermore, at variance with the Lyapunov exponents, it shows jumps when 
there are sudden changes on the available phase-space. 
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One of the most characteristic feature which is emerging when dealing with nonlinear 
systems is the appearance of chaotic motion. There has been a considerable amount of work 
to establish what are the conditions for a nonlinear system (dissipative or conservative) to 
display chaotic dynamics and what are some suitable quantities to characterize it. The 
following gives a brief review of the quantities discussed so far in the literature jT|: 

i) The Lyapunov exponent (LE), the mean rate of separation between two adjacent 
trajectories in phase space, is one of the most used measures. A system is chaotic when the 
trajectories diverge exponentially, i.e. the LE is larger than zero. 

ii) The correlation function characterizes the "memory" along one trajectory. It decays 
quickly to zero in the chaotic regime. 

hi) The power spectrum changes from discrete lines to a broad-band noise when chaos 
sets in. 

These quantities can be calculated from the nonlinear equations of motion which describe 
the system. Often, however the equations of motion are not known but there might be some 
experimental determination of the time evolution of a physical quantity. There are in the 
literature some suggestions on how to extract the LE from a time series |2|. However the 
proposed methods are not unique and depend on some working parameters . Other cases of 
interest exist where there is no information at all about the time evolution of some quantity 
and final phase space distributions are only known. For instance in the field of nucleus- 
nucleus collisions there is an active search for a liquid to gas phase transition at excitation 
energies of tens MeV/particle or a phase transition to a quark gluon plasma at much higher 
excitation energies. In these experiments the final momentum distributions (of almost all) 
particles are known. There is absolutely no information about the time evolution of the 
system and there is no way to estimate the LE from data. Clearly, it becomes very important 
to have hold of some clearly defined physical quantity that could define unambiguously the 
occurrence of a phase transition. It is the purpose of this paper to show that there is a 
quantity, the asymptotic distance between trajectories, which could be easily estimated in 
such cases where the quantities (i-iii) cannot be unambiguously determined and which gives 
for instance as much information as the LE alone. In particular we will describe how the 
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distance between trajectories, which defines the LE, saturates because of the finiteness of 
the available phase space. We would like to stress that in the literature particular care 
has been taken to estimate the LE by avoiding such saturation of the trajectories. Here, 
we are going to do the opposite, i.e. study in detail the saturation properties of trajectories 
to demonstrate how many physical informations can be obtained from such studies. We 
will discuss some illustrative examples using one- dimensional maps. Preliminary studies for 
Hamiltonian systems can be found in Q. 
A one-dimensional map can be written as: 

x n +i = f(x n , r) (1) 

where r is a control parameter. Such maps, even though apparently very simple, exhibit 
most of the chaotic features found in more complex systems. 

The appropriate quantity to discuss here is the mean distance between two trajectories 
separated initially by a very small distance do: 

I n 

^ n = ~N ? ^ Xfl ~ X n ^ ^ 

i=0 

and: 

4 ] = r(4°) (3) 
xf = r (4° + (4) 

In order to avoid fluctuations arising from a particular choice of the starting point for the 
iterations we average over N couples of trajectories and choose the initial starting points x l 
randomly from a uniform distribution. After n iterations, for small enough initial separation, 
the approximate n dependence of the distance d n is: 

d n = e nX d = e x d n _i = Ad n ^ (5) 

where A = e A and A is the Lyapunov exponent of the map: 

i=n 

X = lim -Y\ln | f{x n ) \ (6) 

i=l 
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the prime indicates the derivative of the function f(x n ) respect to x n [p]]. The action of the 
maps consists generally of two steps: the stretching, which leads to the exponential regime 
and is characterized quantitatively by the LE and the folding process which keeps the orbit 
bounded. Therefore when n is sufficiently large, the solution eq.(5) is no longer valid. We 
can consider eq.(5) as a first order expansion in d n , and in the hypothesis that: d n < 1 for 
any n , we include a second order correction term in (|^): 

d n+1 = Ad n - Yd 2 n = F(d n ) (7) 

Note the analogies to the derivation of the logistic map ||. 

Let us define the asymptotic value of the mean distance between two trajectories as: 

I n 

4o = lim - V di (8) 

II i=l 

The fixed points of (0) are d\ — and: 

A-l. 

d 2 = d 00 = — - — (9) 

Thus, (0) describes the irreversible approach of the system to equilibrium which correspond 
to the fixed points solution (9). From the stability condition | F'(di) |< 1 [|I[|, we find that 
d\ is a stable point for A < 0, while d 2 = d^ is a stable point for < A < ln(3). We also 
stress these particular and interesting cases: 

a) F'ldoo) = gives A = ln(2) which is a superstable point of the map F(d n ). Notice that 
such value of the LE is obtained in the triangular and logistic maps for a control parameter 
where the maps become ergodic. Thus ergodicity of the maps is equivalent to a superstable 
point of our proposed application describing the evolution of the distance between trajectories. 
Also, di = becomes a superstable point for A — > — oo. 

b^^cioo) = — 1 gives A = ln(3) which is the value for which the map F(d n ) has a 
pitchfork bifurcation. In this case the values of doo at bifurcation are given by the conditions 
d\ = F(d 2 ) and d 2 = -F(c?i). This and larger value cases of the LE are outside the purpose 
of this paper and will be discussed in a future publication . 

The actual value of T can be easily obtained inverting (H): 
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The entire evolution of the distance between trajectories is given by equation ([F|). It 
contains three characteristic quantities, A, doo and T, but only two of them are indepen- 
dent because of the relation (PH[). To better grasp the meaning of the physical quantities 
introduced above let us consider the (unbound) map: x n+ i = 2x n . It easy to show that 
for this map it is A = ln2, and T = 0. If we impose the condition on the map to have 
modulus 1 (Bernoulli shift), the LE remains the same while T = l/d^ = 3. Thus the LE is 
only sensitive to the stretching mechanism while T is sensitive to the folding and stretching 
mechanisms, eq. QTOf) : its knowledge allows us to distinguish between maps which have the 
same LE. 

In fig.(l) we plot d n versus n as obtained numerically for the logistic, triangular and sin 
maps, for three different initial values of do 0, full lines. We observe in all cases that, after 
a fast increase, the distance between trajectories saturates. The value of saturation is d^ as 
defined in (8) independent on the initial relative distance d^. Inserting the values of A and 
d^, as obtained from eqs.(6) and (8), in eqs. (7) and (10), gives the dot points displayed in 
fig.(l). The agreement to the numerical results is extremely good for all cases and supporting 
our hypothesis. We conclude that in order to have an overall description for the evolution 
of the mean distance between trajectories, beyond the exponential regime cf. eq.(5), we need 
two parameters, A and d^, cf.eq.(7). 

At the ergodic point, corresponding to fully developed chaos, we can calculate analytically 
the LE and d^ as mean values over the phase space using the corresponding distribution 
function |]. We obtain: 

A e = f 1 p( x )ln | 4 - 8x |= ln2 (11) 
Jo 

rl 4 

dooe= p(x)p(y) \ x - y \= -~ (12) 
Jo vr z 

for the logistic map. For the triangular map we have A e = ln2, and d^ = 1/3, the subscript 
e stands for ergodic. These values are in perfect agreement with the ones used in fig. (la), 
and fig.(lc). 



For other values of the control parameter it is not possible to calculate analytically A 
and doo and we have performed some numerical calculations displayed in fig. (2). Figs. (2a) 
and (2c) give the LE and d^ respectively as obtained from eqs.(6) and (8) for the logistic 
map. The right column refers to the results of the triangular map. Recall that the LE for 
the triangular map is given by A = ln(2r). These are some features of particular interest: 

i) at particular values of the control parameter, d^ has jumps which indicate a change in 
the dynamical behavior. Notice, for instance, the jump near r=0.7 for the triangular map. 
This jump, not observed in the LE, is due to the sudden increase of the available phase space, 
because of band splitting bifurcation, see also fig. (2b); 

ii) similarly to the LE, at the transition to regular windows d^ drops to zero. 

Let us study this last point in more detail for the logistic map. Near the transition point 
from order to chaos, the LE behaves like 

Aoc {r-r^f ; f3 = — (13) 

Here is the accumulation value of control parameter for the double period bifurcation 
cascade. This relation can be easily obtained by means of the Renormalization Group Theory 
(RGT) Within the same framework we can easily derive d^ and obtain |J: 

doo oc(r- roe)^ 2 ; v/2 = 1 ^- (14) 

Ino 

In the inset of fig. (2c), we plot the numerical result for d^ vs. r near the critical point. The 
full line gives the power law, eq.(14), in very good agreement to the numerical values. We 
remark that the critical exponent for doo depends on both Feigenbaum constants a and 8, 
while the LE depends on 5 only. We have used the symbol v in analogy to the treatment of 
second order phase transition. Infact, the counterpart of the distance between trajectories 
for maps is given by the variance in momentum space for Hamiltonian systems Near a 
second order phase transition such a variance is proportional to the inverse of the correlation 
length which depends on (T — T c )~ u , where T c is the critical temperature JJ. 

To make our analogy to phase transitions more complete recall that at the critical point 
the correlation function, defined as 0: 

6 




1 ^ 



(15) 



-m 



where ^ = / l (x ) - x ai) and x, 



av 



lim 



■oo 



Sr=i x «> decays with a power law in m: 



C(m, f r J oc m 



(16) 



V 



ln2 



which depend now only on a. 



All above three critical exponents can be related as follows: 



v — (3rj 



(17) 



which could be considered valid for the transition to chaos through double period bifurcation. 
This behaviour is similar for all the maps having period doubling chaotic bifurcation, the 
critical exponents depending on the values of corresponding a and 5 constants. 



We see that this parameter is very large especially near the critical point for the transition 
from order to chaos, and then it becomes almost constant. The value of the constant is -r—, 
eq.(12). Thus the LE is to a good approximation proportional to especially near the 
ergodic point. This fact is especially important for equilibrated physical systems since one 
can deduce the properties of the LE from the variances which are given from thermodynamics 



In conclusion, we have shown that a useful quantity to characterize the occurrence of 
chaos is g^, the value at which the distance between two trajectories saturates in phase 
space. It is complementary to the Lyapunov exponent giving informations about the global 
features of phase space. It signals also when changes in dynamics are having place which are 
not always reflected in the behavior of the Lyapunov exponent. Both parameters are needed 
in order to obtain a correct description of the time evolution for the mean distance between 
trajectories. In analogy to phase transitions, the LE is the order parameter while d^ is the 
inverse of the correlation length. In some experiments both quantities can be measured, 
often, however, like in nucleus nucleus or cluster cluster collisions, it is not possible to follow 
the time evolution of the system and asymptotic quantities, like g^, can only be detected. 
The generalized application for the distance between trajectories proposed in this paper, 



For completeness in fig. (3) we plot the parameter V vs. r for three different maps ]7|. 
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eq.(7), gives a link between the initial exponential expansion and the final equilibrium stage. 
We have shown that ergodicity is simply a superstable point of our proposed eq.(7) and 
naturally gives a Lyapunov exponent A = ln2. 
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FIGURE CAPTIONS 

Fig. 1 Distance between trajectories vs. iteration for different maps, (a) Logistic map with 
the control parameter r=4; (b) the same as (a) but for r=3.771; (c) triangular map 
at r=l;(d) sine map at r=0.73. The different curves correspond to different starting 
values of do- The dots are obtained from eq.(7). 

Fig. 2 Lyapunov exponents and asymptotic distances between trajectories vs. control pa- 
rameter r,for the logistic map (left column) and triangular map (right column). In- 
set (2c): vs. r near the critical point for order to chaos transition, for the logistic 
map . The full lines give the power law dependence of doo as discussed in the text. 

Fig. 3 T vs. r for the logistic (a), sine(b) and triangular map(c). The dashed line gives the 
values of Y when A = /n(2) for the three maps respectively. 
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